Statistical mechanics of homogeneous partly pinned fluid systems 
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The homogeneous partly pinned fluid systems are simple models of a fluid confined in a disordered 
porous matrix obtained by arresting randomly chosen particles in a one-component bulk fluid or 
one of the two components of a binary mixture. In this paper, their configurational properties are 
investigated. It is shown that a peculiar complementarity exists between the mobile and immobile 
phases, which originates from the fact that the solid is prepared in presence of and in equilibrium with 
the adsorbed fluid. Simple identities follow, which connect different types of configurational averages, 
either relative to the fluid-matrix system or to the bulk fluid from which it is prepared. Crucial 
simplifications result for the computation of important structural quantities, both in computer 
simulations and in theoretical approaches. Finally, possible applications of the model in the field of 
dynamics in confinement or in strongly asymmetric mixtures are suggested. 



I. INTRODUCTION 

The physics of fluids under nanoscale confinement is a 
topic of great interest, from a fundamental, applied and 
interdisciplinary point of view at once (see [l|-|3| for re- 
views). This is however a rather difficult one for theorists. 
Indeed, since many porous solids are structurally disor- 
dered, one is usually faced with a complex interplay of 
finite size effects, adsorption phenomena, topology, and 
randomness, which is not easily captured with analyti- 
cally or numerically tractable models. 

In the past few years, a widely used class of mod- 
els for theoretical and computational studies of fluids 
adsorbed in disordered porous media has been one in 
which the fluid molecules evolve in a statistically ho- 
mogeneous random array of particles frozen in a dis- 
ordered configuration sampled from a prescribed prob- 
ability distribution. Recent works include investiga- 
tions of the non-equilibrium phenomena in the adsorp- 
tion/desorption processes 4|-M| , of the critical behavior 
at fluid-fluid transitions [§-fl2| , and of the glassy dynam- 
ics in confinement |13l - [24j . 

The first elaborate theoretical treatment of this type of 
systems has been derived by Madden and Glandt [25| for 
the model of the so-called "quenched-annealed" (QA) bi- 
nary mixture. In this model, the probability distribution 
of the porous matrix is chosen as the equilibrium distribu- 
tion of some simple fluid system, so that the solid samples 
can be thought of as the results of instantaneous ther- 
mal quenches of this original equilibrium system, hence 
the denomination "quenched" for the matrix component. 
Then the fluid (the "annealed" component) equilibrates 
in the random potential energy landscape created by the 
frozen matrix particles. Thanks to the property of sta- 
tistical homogeneity of the solid, Madden and Glandt 
have demonstrated that the QA mixture can be studied 
with great ease via simple and direct extensions of stan- 
dard liquid-state theoretical methods. In their pioneering 
and subsequent works [25l - l28j . they used diagrammatic 
techniques to investigate the distribution functions and 
thermodynamic properties of the system. Later, Stell 
and coworkers rederived, extended, and, in some cases, 



corrected these results with the help of the replica trick 
[2914321 ]. Interesting developments in the framework of 
the density functional theory have been reported as well 

Following this breakthrough, other prescriptions have 
been put forward for the preparation of the disordered 
porous matrix, resulting in a variety of models. In a de- 
pleted system [37| [H[ , the matrix is obtained by quench- 
ing configurations of a one-component fluid and by re- 
moving at random a fraction of its particles. Strate- 
gies for a correlated depletion step have recently been 
suggested [H, Eo|. In a templated system [il| - |44 ]. the 
porous medium is produced by freezing configurations of 
a binary mixture and by removing one of its components, 
which acts as a template for the remaining one. In both 
models, the probability distribution of the solid can be 
explicitly and easily related to the equilibrium distribu- 
tion of an underlying simple fluid system, as in the QA 
mixture. This needs not be the case. Indeed, for many 
aspects of the formalism, it has been shown that it is 
enough to know the structure of the solid, typically at the 
pair level, instead of its full statistics [13, Ea]- This has al- 
lowed theoretical investigations of fluids confined in real- 
istic porous media generated through out-of-equilibrium 
processes, such as the diffusion-limited cluster-cluster ag- 
gregation (DLCA) process [ir3 |. 

In all the above examples, the preparation of the ma- 
trix and thus its statistical properties are independent of 
the fluid that will be imbibed in it. This looks like the 
most reasonable way to proceed, since this actually re- 
produces the usual experimental situation. However, a 
rather popular model in computer simulation studies of 
confined fluids does not display this feature [TH, [l9|, I47F - 
l52j . This model is the randomly pinned (RP) system, 
which is considered in the present paper with one of 
its natural generalizations. Operationally, it can be de- 
scribed as follows. While exploring the configuration 
space of an equilibrium bulk fluid system, snapshots are 
selected at random times, in which a fraction of the par- 
ticles is chosen randomly and pinned down, thus generat- 
ing samples of rigid disordered matrices. The remaining 
unpinned particles then become the confined fluid phase 
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which evolves under the influence of the same interactions 
as in the original bulk system. Note that a very similar 
setup can be realized experimentally in two dimensions, 
by squeezing a binary colloidal mixture between two glass 
plates [53l |54|. 

It will be shown in the following that the RP model is 
actually a special case of a templated system. As such, its 
core physics does not display any fundamental difference 
compared to any other type of fluid-matrix model. But, 
in this specific setup, it is also clear that the fluid and 
matrix properties are unusually and irremediably entan- 
gled. In fact, since the porous matrix is prepared in the 
presence of the fluid, there is a very peculiar complemen- 
tarity between the mobile and frozen phases, from which 
we will show that nontrivial configurational properties 
emerge which make the RP model worth special atten- 
tion. Furthermore, among the different particle-based 
models of disordered porous media, the RP system is 
unique in allowing of an interpretation as a limiting case 
of an equilibrated binary mixture, whereby the mass of 
the pinned particles is sent to infinity (for Newtonian 
dynamics) or their free-diffusion coefficient to zero (for 
Brownian dynamics) [47M491 152| . Its study is thus rele- 
vant for the understanding of the dynamics of mixtures 
with a strong dynamical asymmetry, of which it repre- 
sents an asymptote. 

One of the properties that will be established in the 
present paper has already been observed in the above- 
cited computer simulation studies, namely that the pair 
distribution functions of the bulk system on which the 
RP model is based are preserved after the pinning pro- 
cess [H, IH, [5(| ■ Up to now and depending on the 
authors, however, it was never clear whether this finding 
was reflecting an exact property of the system or was just 
a good approximation. The issue is settled here in favor 
of the first interpretation. In addition, it will be shown 
how other distribution functions, which are characteris- 
tic of systems with quenched disorder and usually quite 
difficult to compute efficiently, should be easily accessible 
in the very special case of the RP model. 

The paper is organized as follows. In Sec. II, the RP 
model and a useful generalization are described with their 
defining probability distributions. In Sec. Ill, configura- 
tional identities are derived, which relate different types 
of averages. They are the main results of this work, whose 
possible applications in computer simulation studies are 
discussed. Their consequences for the structure at the 
pair level and their implications for integral equation the- 
ories are considered in Sec. IV, while a dynamical point 
of view is developed in Sec. V. Section VI is devoted to 
concluding remarks. 



II. MODELS AND PROBABILITY 
DISTRIBUTIONS 

In this section, the RP model is defined and its es- 
sential statistical properties are derived. To achieve this 



goal, we proceed in two steps. First, a formal link is es- 
tablished between the RP and templated systems, which 
leads to a natural generalization of the problem. Second, 
the relevant probability distributions are provided in this 
extended framework. 

The simplest version of the RP model, which is suffi- 
cient to illustrate the main properties of this type of sys- 
tems, derives from a one-component bulk fluid through a 
random pinning process, as sketched in Fig. Q] (55[. More 
specifically, for each configuration of the fluid which oc- 
curs with a probability distribution corresponding to a 
prescribed statistical ensemble, a fraction x of its parti- 
cles is randomly chosen and pinned down to form a disor- 
dered porous matrix, while the complementary \ — x frac- 
tion remains mobile and provides the confined fluid. Dif- 
ferent ensembles can be considered, which are all equiv- 
alent in the thermodynamic limit l32f . In existing com- 
puter simulation studies [f8|, [l9|, |47I452| . both the fluid 
and the matrix have always been treated in the canonical 
ensemble, i.e., with no fluctuations in their particle num- 
bers and thus in the pinning fraction. Here, we find more 
convenient to resort to a grand-canonical description in 
which, in particular, the pinning fraction is allowed to 
fluctuate around its mean value. 

Accordingly, we start with a one-component bulk fluid 
in a volume V at temperature T (as usual, we define 
j5 = \/k^T) and activity z. From well-known statisti- 
cal mechanics [56}, the configurational probability den- 
sity of finding this system with N particles located at 
(ri,r 2 , ...,rjv) 
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where Vbuik(-V, ? N ) denotes the potential energy of the 
fluid. The partition function Sbuik is simply 



^bulk — Tr 



(2) 



where the shorthand notation for the trace operator 
Tr '-- = EN=oJdr N -- - is used. 

During the pinning process, for every configuration of 
the fluid, each particle in the system might be pinned 
down with probability x or left mobile with probability 
I — x. It results that the joint probability density of 
generating a matrix with N m immobile particles located 
at (qi,q2, ■ • ■ ,qw m ) = q"" 1 while N t unpinned particles 
are located at (sj.,S2,...,Sjv t ) = s Nt a t the time of the 
pinning process, with N m + N t — N , is 



Pmt(^ m ,q W »,V t ,s iVt ) = 



\N t e -l3V b ^(N m +N t ,q" 



SbulkV m W t ! 



(3) 



which is basically Eq. (|T|) modified by a combinatorial 
factor due to the random pinning process. For clarity, 
different notations are used for the locations of the pinned 
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(a) 



(b) 




FIG. 1. Simple schematics of the fluid-matrix models studied in this work. In all panels, the immobile particles are represented 
in black, and the mobile ones in white, with arrows symbolizing their movement, (a) Randomly pinned system. In a precursor 
one-component fluid (left panel), a predefined fraction of the particles is chosen at random in an instantaneous configuration 
and pinned down to form the matrix (middle panel, the pinning fraction is x — 1/3), in which the remaining particles continue 
their motion (right panel), (b) Partly pinned system. Starting with an equilibrated binary mixture (left panel), one component 
is pinned down in an instantaneous configuration to form the matrix (middle panel), in which the other component continues 
its movement (right panel). 



and unpinned particles, and, in anticipation of the com- 
ing interpretation, the subscripts m and t for "matrix" 
and "template" , respectively, have been introduced. The 
probability density of the matrix configurations immedi- 
ately follows by tracing out the unpinned particles, 

V m (N m , ) = Tr t V mt (N m , , Nt , S Nt ) . (4) 

This is exactly what would be obtained for a depleted 
system (37l. l38j. This is not surprising, since it is rather 
intuitive that a RP system is equivalent to a depleted 
system in which the depleted particles are reinjected as 
the confined fluid. 

A similar equivalence can be established with a tem- 
plated system (4lTj44l |. Indeed, thanks to the binomial 
identity, Sbuik can be rewritten as 

Sbulk = Tr m Tr t 

,JV»+N, T JV m /l _ T Wt p -/3Vbuik(JV m +Af t ,q JV '«,s Jv t) 



Pmt^m, q^ 1 ", N t , s Nt ) is then easily recognized as the 
grand-canonical configurational probability density of an 
ideal binary mixture, with activities z m = zx and z t = 
z(l — x) for the components m and t, respectively, and 
V m (N m , q^-) as the probability distribution of the tem- 
plated matrices that can be generated from this ideal 
mixture, components m and t being the matrix and tem- 
plate components, respectively. Note that Vbuik appears 
unchanged, because the potential energy of an ideal bi- 
nary mixture is, by definition, independent of its compo- 
sition and thus equal to that of any one of its components 
alone. Therefore, a RP system with a pinning fraction x 
is equivalent to a templated system based on an ideal bi- 



nary mixture with a matrix number fraction x, in which 
the template particles are reinjected as the mobile phase. 

An immediate generalization of the problem follows 
from this equivalence, by relaxing the restriction to ideal 
mixtures (see Fig. [1]). We call the corresponding ex- 
tended class of models the partly pinned (PP) systems, 
which might be studied generically and only specialized 
afterwards to deal with the RP model. Thus, instead of 
Eqs. ([3]) and (|5|), we shall consider the generic expressions 
for templated systems, i.e., 

P m t{N m ,q N ™,N t ,s Nt ) = 

_ JV m _ N t -p [y mm (JV m ,q N » )+V mt (JV m ,q"» ,N t ,s N t ) + Vtt (jv t :S «t )1 
z m z t e _ 

SbuikAUWt! 

(6) 

for the grand-canonical configurational probability dis- 
tribution of the bulk matrix-template system and for the 
joint probability density of generating a matrix with N m 
particles located at q Wm while Nt template particles are 
located at s^*, with 

Sbuik = Tr m Tr t 

?N m _ Nt -P [y mm (N m .q™'" )+V mt (A^q™'" ,Nt ,s"t )+Vtt (N t ,s N * )1 
2 m z t e 

NJ.N t \ 

(7) 

for the normalizing partition sum. Equation Q still gen- 
erates the matrix probability distribution. In the above 
equations, z m and z t are the activities of the matrix and 
template components in the bulk binary mixture, and 
Kn m (N m , q^ ) , V mt (N m , q^ ,N t ,a N *) and V tt (N t ,s N <) 
are the matrix-matrix, matrix-template and template- 
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template contributions to the potential energy of the mix- 
ture, respectively. 

So far, only probabilities relative to the matrix and its 
preparation process have been considered. We now turn 
to the confined fluid statistics. As mentioned above, the 
confined fluid is assumed to be in grand-canonical equi- 
librium in the presence of the porous matrix and, by con- 
struction of the model, to inherit its properties from the 
template. Thus, its activity Zf (the subscript f is used for 
properties of the confined fluid) is taken equal to the tem- 
plate activity z t and, for a particular realization of the 
matrix with A m particles located at q Nia , its potential 
energy is chosen as V mt (N m , q N ™,N h r N <) + V tt (N h r N <), 
when Nf fluid particles located at (n, ra, . . . , rjv f ) = r Nl 
are present in the system. It follows that the probability 
density of such a configuration, which depends paramet- 
rically on the matrix configuration (A^q^" 1 ), is 



(8) 



Pf(JV f ,r w *|JV ni ,q- 

z Ar fe -/3[y mt (Ar m ,q JV »,iV t ,r N f) + ^t(A f f,r N f)] 
S f (A m ,q^)A f ! 

with the confined fluid partition function 



E f (A m ,q w ») = 

N t -P \ V m t ( , q Nm , N t , r N f ) +V tt ( N t , r N i )] 

7S ' (9 » 

Combining Eqs. (J6j) and ((SJ, a very simple equality can 
be derived, 

Vmt (N m , q^" , JV t , )V { (Nf, v N < | N m , q"- ) = 

7> m t(A m , q N -,N f , r Nt )V { (N t , b n * \N m , q*»), (10) 

which reflects the very peculiar symmetries of the system 
under study and will play a crucial role in the following. 

Armed with these results, we might now investigate the 
configurational properties of the RP and PP systems. 



III. CONFIGURATIONAL IDENTITIES 

Dealing with quenched-disordered fluid-matrix sys- 
tems, two types of configurational averages have to be 
considered when computing their properties [III [27j • The 
first one is the ordinary thermal average denoted by (• • • ), 
taken for a given realization (A^q-^™) of the matrix 
with the probability density Vf(Nf, r Nl | A m , q^'"). The 
second one is the disorder average over the matrix real- 
izations, denoted by 7TT , to be taken with the probabil- 
ity density V m (N m ,q Nxa ) after the thermal average. We 
shall not dwell on subtleties of the physics of quenched- 
disordered systems, but it might be worth recalling that 
this average over disorder is equivalent to an average over 
macroscopic subparts of a macroscopic system for addi- 
tive quantities, and to a volume average over a macro- 
scopic sample for locally defined quantities 57 1. 



These two types of averages might be combined in 
many different ways. For instance, computations of free 
energy differences typically involve expressions of the 
form ln{A). Here, we shall concentrate on two specific 
examples for which simple results can be derived. They 
correspond to double averages (A) and to products of 
the form (A){B). Such quantities are often combined 
to generate correlation functions characterizing the two 
physically distinct types of fluctuations present in dis- 
ordered systems. Thus, the typical thermal fluctuations 
are qu a ntified by so-called connected averages of the form 
(AB) — (A)(B), while the disorder-induced fluctuations 
of thermal quantities are measured b y so-call ed discon- 
nected averages defined as (A)(B) — (A) (B). The dis- 
tinction between these two types of correlation functions 
is an essential feature of the physics of disordered sys- 
tems, whose significance is, for instance, stressed by the 
fact that the thermodynamic susceptibilities, such as the 
isothermal compressibility [28l |32| , are always expressed 
as connected averages. 

In order to shorten the equations in this rather for- 
mal section, a condensed vector notation will be used. 
Thus, we define m = (N m ,q N ™), t = (N t ,s Nt ), and 
f = (A f , r Nt ). With these definitions, Eq. (|TU)) now reads 



P mt (m,t)P f (f|m) = P mt (m,f)P f (t|m). 



(11) 



Generically, a configurational variable for a quenched- 
disordered fluid-matrix system is a function A(x; m, f ) of 
the matrix and fluid particle numbers and coordinates, 
possibly with other variables (such as space variables 
when dealing with n-particle densities) collectively de- 
noted by x. Its thermal average for a given matrix real- 
ization is defined as 

(A(x; m, f)> = Tr f Vt(f |m)A(x; m, f), (12) 

so that its double average after tracing out the matrix 
variables reads 



(A(x;m,f)) =Tr m P m (m)Tr f P f (f|m)^(x;m,f). (13) 

Specializing to the PP systems and introducing an ex- 
plicit reference to the original matrix-template mixture, 
this can be rewritten as 



(A(x; m, f)) - Tr m Tr t Tr f P mt (m, t)V { {f |m)A(x; m, f), 
which, thanks to Eq. (|11[) . is transformed into 



(A(x;m,f)) = Tr m Tr t Tr f P mt (m,f)P f (t|m)A(x;m,f), 

(15) 

eventually leading to 



(A(x;m,f)> =Tr m Tr f P mt (m,f)A(x;m,f), (16) 



since Trt 'Pf(tlm) = 1 appears factorized in Eq. (|15|) . 

The right-hand side of Eq. (|16[) is immediately recog- 
nized as a thermal average for the bulk matrix-template 
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mixture, which will be denoted by (• • • )bulk- So, the first 
sought-for configurational identity reads 

(A(x;m,f)) = (A(x; m, f)) bulk , (17) 

i.e., the double-averaged configurational properties of a 
PP system coincide with the corresponding quantities in 
the bulk fluid on which it is based 58] . This in particular 
applies to the n-particle densities and distribution func- 
tions and ex plai ns p ast observations in computer simula- 
tion studies [H SI |H HI! ■ 

There are different possible ways of taking advantage 
of this result in computer simulation studies. An immedi- 
ate idea is that one can completely avoid the computation 
of double averages for the PP fluid-matrix systems. In- 
stead, the corresponding calculations can be performed 
on the bulk fluid, which in any case has to be simulated 
in order to generate the porous samples and for which 
the problem is conceptually simpler (only one type of 
average is required and it is an ordinary thermal aver- 
age) and the sampling of the configuration space is often 
more efficient, thanks to faster dynamics [la. l47l l48l l50l- 
[52( | . Another option is to use the computation of dou- 
ble averages as a means to calibrate the parameters of 
the simulation study. Indeed, the above identity is not 
expected to hold on a sample-by-sample basis. It be- 
comes valid only after the disorder average is performed. 
Therefore, one can use the comparison between double 
averages and high quality data for the bulk as a guide to 
estimate the minimal number of matrix realizations that 
is required in order to achieve a satisfactory convergence 
of the disorder averaging procedure. 

In the case of two configurational variables A(x; m, f ) 
and B(y; m, f), the typical value of the product of their 
thermal averages is defined as 

(A(x;m, f))(B(y;m, f')> = Tr m P m (m) 
Tr £ P f (f |m)A(x; m, f ) Tr' { V { (f'\in)B(y; m, f ), (18) 

where Tr' { simply represents the trace over primed fluid 
variables. Through Eq. (QJ, this might be rewritten for 
a PP system as 

{Kx; m, f))(B(y;m, f')> = Tr m Tr t Tr f Tr^ V m t (m, t) 
V f (f |m)A(x; m, f )7> f (f |m)B(y; m, f ). (19) 

Thanks to Eq. (fTTj) , one can exchange the variables t and 
f between 7> m t(m,t) and Vi(f'\m) (one could choose to 
exchange t and f instead, but this would simply amount 
to an exchange of A and B which play symmetric roles in 
the original problem), then perform the trace Tr t which 
reduces to a normalization condition as above. It remains 

(A(x; m, f))(B(y;m, f')> = Tr m Tr f TrJ 

V m t (m, f)V f (f |m)A(x; m, f )B(y; m, f), (20) 

whose meaning is illuminated by a simple change of 
dummy variables leading to 

(^(x; m, f))( J B(y;m, f')> = Tr m Tr t Tr f 

P mt (m, t)P f (f |m)A(x; m, f)B(y; m, t). (21) 



The right-hand side of the latter equation takes the 
form of a double average, but now with a modified dis- 
order average involving the matrix-template probability 
distribution V m t(m, t), which will be denoted by 7TT ' in 
the following. Such a modification is very natural if one 
includes the template variables ( N t , s Nt ) in the set of 
configurational parameters describing the system and ac- 
cordingly deals with configurational variables of the form 
A(x;m, t,f). Strictly speaking, this extension is not re- 
quired by the physics of the system, which can be dis- 
cussed uniquely in terms of fluid and matrix parameters 
[27| . in which case, as it should be, the modified disor- 
der averaging procedure does not change anything, as 
attested by identities such as 

(A( X ; m,f))' = (A(x;m,f)). (22) 

However, the consideration of template degrees of free- 
dom offers a very pictorial way to capture how the 
matrix-template and template-template correlations im- 
printed in the porous solid during its preparation are 
transferred to the fluid. It can also lead to important 
technical simplifications in formal developments. One 
example based on the interpretation of Eq. (f2Tj) will fol- 
low. Another one is provided by the integral equation 
theory of the generic templated systems [Tll-l44l|. where 
the introduction of total and direct correlation functions 
involving the template leads to equations that preserve 
standard diagrammatic prescriptions (no nodal points in 
direct correlation functions) , at variance with Madden's 
more compact formalism [27]. 

Coming back to Eq. (l2"Tj) , the second sought-for con- 
figurational identity thus reads 

(A(x; m,f))(B(y;m,f) = (A(x; m, f )B{y; m, t))', 

where, since the matrix and template variables are static 
quantities, £?(y;m,t) can be placed inside or outside 
the thermal average. So, in a PP system, the disorder- 
averaged product of two thermal averages coincides with 
the modified double-averaged product of a fluid-matrix 
configurational variable with a template-matrix function. 
Interestingly, this result can be formulated in terms of 
non-equilibrium quantities as well. Indeed, since the fluid 
and the template components are essentially identical, 
the template coordinates at the time of preparation of 
the matrix can be considered as initial conditions for an 
ulterior fluid dynamics inside the matrix. So, the right- 
hand side of Eq. (|2U|) can be interpreted as the correlation 
of a thermal average with a function of these initial con- 
ditions. 

The availability of Eq. (|23[) represents a major simplifi- 
cation for computational studies of PP systems compared 
to other particle-based models of disordered porous me- 
dia. Indeed, the modified double average on its right- 
hand side should be no more difficult to compute than 
the standard double average on fluid and matrix con- 
figurations met in all this family of systems. The only 
additional price to pay is to keep track of the positions 
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occupied by the template particles at the moment of the 
matrix production. Thus, through Eq. (|2"5)l , quantities 
of the form (A){B), which are required to compute the 
physically important connected correlations, are made 
numerically accessible via a simple and controlled pro- 
cedure. 

The relevance of this result is best illustrated through 
an example, which is the computation of the two-point 
correlation function 

V(r, r') = (pUfofMpW^f')), (24) 

where pp (r; f ) is the microscopic fluid density operator 

N, 

4 1) (r;f)=^^(r-r l ). (25) 
i=l 

For a generic fluid-matrix system, one typically has to 
calculate (pP(r;f)) on a grid for each matrix realiza- 
tion and to compute the autocorrelations as sums over 
the grid points, before disorder averaging [59| . Unfor- 
tunately, the outcome of this procedure has been found 
to be marred by artifacts related to the finite size of the 
grid cells. Analogous difficulties are present in recipro- 
cal space calculations of the structure factor associated 
with -0(r, r') [goj]. In the case of a PP system, however, 
Eq. (T2"5)l can be used to obtain 

V(r, r') = (^(r^pfVst))', (26) 

where p^'(r;t) is the template analogue of pfjpfcf). 
Further simplification occurs by noting that the prod- 
uct p]p (r; f )p[ 1 ' (r'; t) is simply the two-particle fluid- 

template density operator p ft ; (r, r'; f , t), so that even- 
tually 

^(r,r') = (pf (r,r';f,t))''. (27) 

The computation of ip(r, r') is thus turned into the rather 
innocuous problem of computing the pair correlation 
function between the fluid particles and the template 
sites, which does not require the use of a grid and thus 
is obviously free of grid artifacts. It is interesting to note 
that an analogous possibility to replace the direct compu- 
tation of a disconnected correlation function by the com- 
putation of an equivalent pair correlation function also 
exists for ideal gases in disordered environments, thanks 
to the fact that the connected pair correlation functions 
identically vanish in these systems. The case has been 
investigated on one specific example in Ref. [59| |. where 
it is reported that the latter strategy yields better results 
than the former. This clearly lends support to the sug- 
gestion that the availability of Eq. (|2"5|) and the resulting 
simplifications for PP systems should facilitate the accu- 
mulation of high quality computer simulation data. 

Finally, the above equations, formulated in the frame- 
work of the PP systems, can be easily adapted to the case 



of the one-component-fluid-based RP models. In particu- 
lar, using Eq. (fTT|). the double-averaged quantities in the 
RP system can be related to corresponding thermal aver- 
ages in the one-component bulk fluid. One only needs to 
revert the mapping to an ideal template-matrix mixture 
described in the previous section. Unfortunately, the sim- 
plicity of the above identities is lost in the process, since 
the dependence of the resulting expressions on the pin- 
ning fraction is found to change with the configurational 
variable under consideration. For this reason, we do not 
report any specific results for the RP model here. Some 
important ones will appear in the next section. 



IV. PAIR CORRELATIONS AND INTEGRAL 
EQUATION THEORIES 

The description of the structure of fluids at the pair 
level is a central issue of liquid state theory. In this sec- 
tion, we examine how the above identities impact on this 
problem for the PP and RP fluid-matrix models. 

Since we are dealing with special cases of templated 
systems, the natural framework for the present study 
is the theory developed by Van Tassel and coworkers 
[4"U - |44j ] . For a generic templated system with no spe- 
cial symmetries, they introduced eight total correlation 
functions in order to fully describe the structure at the 
pair level. Six of them are the standard functions re- 
flecting the two-body correlations between the differ- 
ent types of particles involved in the model. They are 
the matrix-matrix, matrix-template, template-template, 
matrix-fluid, template-fluid, and fluid-fluid total correla- 
tion functions, denoted by h min (r), h mt (r), h tt (r), h mi (r), 
ht{(r), and /iff(r), respectively. They can be expressed in 
terms of double-averaged two-particle density operators 
Plj ; (r,r'), i,j = m,t,f (for notational convenience, the 
dependence on the particle numbers and configurations 
is omitted here and in the following), as 

hij (\r-r'\) = W^lzfiEi , (28) 
PiPj 



where pi = (fi\ (r)) is the number density of species i 

[see Eq. (T2"5|) for the expression of (r)]. Note that a 
modified (primed) disorder average is used, because it is 
required for the proper definition of quantities pertaining 
to template particles. When only matrix and/or fluid 
particles are involved, it can safely be replaced by the 
usual matrix average on the basis of Eq. (|22|) . Through 
these averages, the statistical isotropy and homogeneity 
of the system are restored, hence the dependence of all 
pair correlation functions on the scalar |r — r'| only. 

The two remaining functions are specific to quenched- 
disordered systems and give an account of the correla- 
tions between the one-body fluid densities at two differ- 
ent points induced by the disorder. They are the blocking 
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or disconnected total correlation function 



and the connected total correlation function 



(29) 



hc Qr-^=^'»-f^^, (30) 

Pi 

hence the discussion of the function tp(r, r') in the pre- 
vious section. Note that /tff(r), /i b (r) and h c (r) are not 
independent, since they obey ha{r) = hb(r) + h c (r). 

Using diagrammatic techniques or the replica trick, 
Ornstein-Zernike (OZ) equations are easily derived which 
connect the above total correlation functions to the corre- 
sponding set of direct correlation functions, (kj(r), Cb(r) 
and c c (r). They read (for convenience, the r dependence 
is omitted) 

+ Pm Cmm ® ^ mm + PtCmt ® h mt , (31a) 
Cmt ~l~ PmCmm ® h mt + ptCmt ® /itt, (31b) 

fott = C tt + PmCmt ® /l mt + PtCtt <8> /j-tt, (31c) 
® ft-mf + jOtCmt ® ^tf + PfCmf <8> /l c , 

(31d) 

/itf = C t f + PmCmt ® h m { + p t Ctt <8> ft-tf + PfCtf <8> /lc, 

(31e) 

ha = c f f + /O m c mf ® /i mf + p t c t f ® /i t f+ 

PfCff ® /iff - pfc b ® /ib, (31f) 
/ib = c b + p m c mf ® h mi + p t c t f <g> /itf+ 

Pic c ®hb + p&b ® /i c , (31g) 
h c =c c + p{C c <g> /i c , (31h) 

where ® represents a convolution in real space. As for the 
total correlation functions, one gets Cff(r) = c b (r) + c c (r). 
The first three equations are simply the OZ equations for 
the bulk binary matrix-template mixture from which the 
porous medium is formed, while the others relate to the 
structure of the fluid adsorbed in the random solid. 

Specialization to the present systems of interest is 
achieved by working out the consequences of the iden- 
tities (fTT)) and (|23|) on the one- and two-body densities. 
For a generic PP system, one obtains 



Pt = Pi, 
h m %{r) = h m {(r), 
fht(f) = hff(r), 
h t {(r) = h h (r), 



(32a) 
(32b) 
(32c) 
(32d) 



and equivalently for the direct correlation functions, so 
that the OZ equations reduce to 



^mm 


f-mm 


+ PmCmm ® h mm + p{C m { (g) /l mf , 


(33a) 


h m i — 


Cmf " 


f PmCmm ® ft m f + PfCmf ® /iff, 


(33b) 


/iff = 


Cff + 


- PmCmi ® Vf + PfCff ® /iff, 


(33c) 


h h = 


c b + 


PmCmf ® /lmf + PfCff <8> /iff - /9fC c <S 


3 /i c , 
(33d) 


h c = 


Cc + 


PfC c (g) /l c . 


(33e) 



For a RP system with pinning fraction x based on a 
one-component bulk fluid with density p and total and di- 
rect correlation functions h(r) and c(r), well-known prop- 
erties of ideal binary mixtures lead to 



/iff(r) 



Pm = xp, 
Pi = (1 - x)p, 
hmi(r) = hm m (r) 



Hr) 



(34a) 
(34b) 
(34c) 



and analogously for the direct correlation functions, so 
that the OZ equations become 



h = c + pc ® h, 
hh = Cb + pc Cg) h — (1 — x)pc c 
h c = Cc + (1 — x)pc c <g> h c . 



(35a) 
(35b) 
(35c) 



Some features of the above equations are very appeal- 
ing with a view to developing integral-equation theoreti- 
cal approaches. The most obvious one is their mere num- 
ber. Indeed, not surprisingly, a significant reduction in 
the number of unknown functions and OZ equations link- 
ing them ensues from the special symmetries of the PP 
and RP systems. It is particularly drastic in the case of 
the RP system, where only two linearly independent OZ 
equations remain [remember that /iff(r) = /i b (r) + h c (r) 
and Cff(r) = c b (r) + c c (r)]. Since one generically needs 
at least two of the three functions /iff(r), hb(r) and h c (r) 
in order to characterize the pair structure of a fluid in a 
statistically homogeneous disordered environment, this is 
actually the smallest possible number, whose only known 
realization so far was in the case of a Gaussian random 
field [6l| . i.e., of a non-particle-based disordered envi- 
ronment. Thanks to this extreme compactness of the 
formalism, with its minimal number of unknowns and 
coupled equations, one can reasonably expect that the 
development and the implementation of advanced inte- 
gral equation theories for the present class of systems 
will be made much easier. 

The situation is actually even more favorable than 
that. Indeed, from Eqs. (|33[) and (|35p . one can imme- 
diately see that the problem of computing the pair cor- 
relations between the different types of particles present 
in the system, which amounts to a calculation on the 
bulk fluid from which the confined fluid-matrix system 
is prepared, is completely separated from the one of ob- 
taining the two contributions to /iff(r) which split under 
the effect of randomness. Therefore, since well estab- 
lished integral equation schemes exist which allow one to 
successfully tackle the first problem for a variety of inter- 
action potentials [56, 62], one can build on them, so that 
the only actual issue for the present systems of interest 
boils down to the computation of hb(r) and h c (r), ha(r) 
being considered as known. 

There arc different possible ways of taking advantage 
of this result. One of them is to note that, because of the 
separation of the bulk- and confinement-specific parts of 
the calculation, one should not feel bound to treat both 
problems with the same type of approximation. There- 
fore, for a given description of the bulk fluid system, one 
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FIG. 2. Structure factors of a hard-sphere randomly pinned 
fluid-matrix system, obtained analytically in the Percus- 
Yevick approximation. Su(q), S c (q), and Sb(q) are, respec- 
tively, the fluid-fluid, the connected and the blocking struc- 
ture factors of the confined system, S(q) is the structure 
factor of the bulk fluid from which the fluid-matrix system 
is prepared. The compacity of the original bulk system is 
4> = (7r/6) pdP = 0.36, with d the diameter of the spheres, and 
the pinning fraction is x — 2/3. 



can freely experiment with closures for hh(r), h c (r), Cb(r) 
and c c (r) . In combination with quantitative comparisons 
with computer simulation data, whose accumulation is 
greatly facilitated by the use of Eqs. (fTT)) and (|2"3"|) (see 
the discussion in the previous section), this should lead 
to insight into this crucial but not so well understood 
aspect of the physics of fluids in random environments, 
with a real prospect that accurate theories can be devel- 
oped, not only for the PP and RP systems, but for fluids 
in quenched disorder in general. 

An alternative strategy is to use the same closure ap- 
proximation for both parts of the calculation and to con- 
centrate on integral equation theories that have analytic 
solutions in the bulk, trying to generate analytic solu- 
tions for the confined system as well. Such solutions are 
extremely scarce in the theory of fluids in disordered en- 
vironments, even within basic approximation schemes. 
They might be useful in a number of different ways. For 
instance, as demonstrated by many advances in the case 
of bulk systems [62[, they might represent precious 
starting points for the development of more elaborate 
approximations. Or they might be used to produce at a 
minimal computational cost the large amounts of struc- 
tural data that are required as input by other types of 
calculations, such as the prediction of dynamical phase 
dia gram s in the framework of the mode-coupling theory 

mm- 

It turns out that closures which fall in the class of 
the so-called Madden-Glandt approximations, i.e., char- 



acterized by an identically vanishing blocking direct cor- 
relation function (29rj3l|. are particularly well suited for 
such an approach. They include the Percus-Yevick (PY) 
and mean-spherical (MSA) approximations, for which, 
precisely, a number of analytic solutions for bulk sys- 
tems are well known (5(| [62[. As an illustration of this 
procedure, we report in Fig. [2] the structure factors of 
a hard-sphere RP system obtained analytically with the 
PY closure. Since Cb(r) = and c c (r) = c(r) in this ap- 
proximation, they are simply expressed in terms of the 
analytically known Fourier transform of the PY direct 
correlation function of the bulk hard-sphere fluid, c(q) 
[f(q) denotes the Fourier transform of f(r)], as 



S(q) = l + ph(q) 



Sti(q) = 1 + Pthff(q) = 
S c (q) = 1 + pth c (q) = 

Sh(q) = Pihb{q) 



1 - pc(q) 
1 - xpc(q) 



1 - pc(q) ' 
1 

1 - (1 -x)pc(qY 
x(l - x)p 2 c{q) 2 



[1 - pc(q)][l - (1 - x)pc(q)] 



(36a) 
(36b) 
(36c) 
, (36d) 



where S(q) is the structure factor of the bulk fluid from 
which the fluid-matrix system is prepared, and Sff(q), 
S c (q), and Sb(q) are, respectively, the fluid- fluid, the con- 
nected and the blocking structure factors of the confined 
system. 

Finally, beside integral equation theories, other liquid- 
state theoretical approaches could either be directly ap- 
plied to the PP and RP systems or at least take advan- 
tage of their simplifying features in some parts of their 
development. For instance, one might turn to pertur- 
bation theory (56| and, in the same spirit as in previ- 
ous work on QA systems [63[ and fluids in aerogels fjfjj ]. 
study athermal PP or RP reference systems decorated 
with attractive fluid-matrix and fluid-fluid interactions. 
One would then benefit from all the above results for the 
description of the reference system. Note in passing that, 
while integral equation theories enforce the complemen- 
tarity between the fluid and the matrix at all state points 
and thus deal with a temperature-dependent structure of 
the porous solid, this would not be the case in such a 
perturbative scheme, where the disordered matrix exclu- 
sively inherits its properties from the reference athermal 
system. 



V. APPLICATIONS IN DYNAMICAL STUDIES 

So far, the discussion has essentially been on how the 
knowledge of the properties of the bulk fluid on which a 
PP system is based leads to considerable simplifications 
in the study of the confined fluid system. Then, one 
might wonder whether, conversely, information on the 
latter might help understanding aspects of the physics of 
the former. In this section, we show that such opportu- 
nities might exist when dealing with dynamics. 
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Operationally, any fluid-matrix system in which the 
porous solid is represented by fixed randomly placed par- 
ticles can be described in dynamical terms as a binary 
mixture of fluid and matrix particles with, respectively, 
finite and infinite masses when the dynamics is Newto- 
nian, or finite and zero free-diffusion coefficients when it 
is Brownian. But, in the case of the PP systems, there is 
more than that. Indeed, it is well known that, in classical 
statistical mechanics, the configurational properties of a 
bulk binary mixture are independent of the finite masses 
or free-diffusion coefficients of its two constituents. Equa- 
tion p7|) shows that this remains true when one mass is 
sent to infinity or one free-diffusion coefficient to zero, 
provided the averaging over the configurations of the now 
immobile matrix particles is performed as a disorder av- 
erage or, equivalently, as a volume average over a macro- 
scopic sample. So, a PP system can effectively be con- 
sidered as an asymptotic case of a binary mixture with 
strong dynamical asymmetry, obtained via a well defined 
limiting procedure [47H49I . [52| . This is a unique situation 
among the different particle-based models of disordered 
porous media, which originates in the fact that the solid 
matrix in a PP system is not prepared independently 
of the adsorbed fluid, but in its presence in equilibrium 
conditions. 

We might now try and transport the above limiting 
procedure into the realm of dynamical theories. In order 
to do so, our starting point will be the description of the 
relaxation of the collective density fluctuations in bulk 
binary mixtures in the framework of the Mori-Zwanzig 
formalism (5(|, first for Newtonian dynamics. 

So, we start with a bulk binary mixture with compo- 
nents labeled f and m. Their particle numbers, number 
fractions, number densities, and masses are denoted by 
Nf and N m , Xf and x m , pf and p m , mf and m m , respec- 
tively. The dynamical variables of interest are the Fourier 
components of the microscopic densities at wave vector 
q and time t, /Of(q, t) and p m (q, t), from which one forms 
the static structure factors, 



Sij(q) 



( Pi (q,0) Pj (-q,0)) 



and density fluctuation autocorrelation functions, 
(Pi(q,t)ft(-q,0)) 



Nt + N n 



(37) 



(38) 



Using standard projection-operator methods [56j, the lat- 
ter functions are shown to obey generalized Langevin 
equations which read in matrix form 

F(q,t)+n 2 (q)F(q,t)+ f drM(q, t - r)F(q, r) = 0, 



(39) 

with initial conditions F(q, 0) = S(q) and F(q, 0) = 0. 
M(q,t) is the matrix of memory kernels and the fre- 
quency matrix ft 2 (q) is given by 



ln% J {q)=q^x l [S-\ 3 {q), 



(40) 



where S -1 (<7) is the matrix inverse of S(q). 

In the limit m m — > +oo, which generates a PP system, 
p m (q, t) is not a dynamical variable anymore, so that 
F m f(q,t)[= Ff m (q,t)] and F mm (q,t) remain equal to the 
corresponding structure factors at all times j64|. As a 
result, the dynamical equation for Ftt(q, t) obtained by 
expanding the matrix products in Eq. (1391) reduces to 



Ff f (q, t) + nUq)Fff(q, t) + n 2 m (q)S m{ ((?)- 



dTMff(q,t-T)F { f(q, T )=0, (41) 



which at infinite time leads to 



n 2 {{ (q)fn(q) + fi? m (9)5 mf (g) +m ff (g) [/ ff (<?) - S fl {q)] = 0, 



where 



fii(q) 



lim lim 

t— too m m — 



Fff(q,t), 



rriff(q) = lim lim Mff{q, t). 

i— >-oo m m — f+oo 



(42) 

(43) 
(44) 



Note that the order of the limits is essential. 

Our goal is to compute mff(q), so it remains to evaluate 
fff(q). If one assumes that the dynamics of the PP system 
is ergodic, it is expected on general grounds that 



fif(q) = lim 



(Mq,*)pf(-q,o)) 



Nf + N n 



(p f (q,i)>(M-q,o)> 



Nf + N B 



XfPfh h (q), (45) 



while Sff(q) = Xf[l+pfhff{q)]. Note that, compared to the 
previous section, additional factors Xf appear due to the 
use of a different normalization for the structure factors. 
By combining Eq. (|4"2"j) with the set of OZ equations (|33p . 
it then results that 



mtf{q) = q 2 



< k B T 



-Pfcb(g)- 



(46) 



For Brownian dynamics with free-diffusion coefficients 
and D° , one gets, instead of Eq. (|3"9"|) . the matrix 
equation [65j 



F(q,t) + fl(q)F(q,t)+ drM(q, t-r)F(q, r) = 0, (47) 
Jo 

where M(q,i) is now the matrix of irreducible memory 
functions and the damping matrix £l(q) is given by 

a ij (q) = q 2 D9x i [S- 1 ] ij (q). (48) 

Repeating the above steps with — > 0, one obtains 

mtt(q) =q 2 D° P fd h {q). (49) 



Equations (|4"6")) and (l4"9l represent bridges between the 
statics of the PP systems and the dynamics of the fully 



10 



annealed mixtures with a strong dynamical asymmetry. 
We suggest that they could be usefully incorporated into 
the development of phenomenological expressions for the 
memory kernel Mff(q, t) reflecting the separation of time 
scales between the short-time relaxation due to the mo- 
tion of the fast particles constrained by the nearly im- 
mobile slow ones and the long-time relaxation associated 
with the ultimate decorrelation of the positions of the 
latter. For instance, since a number of observables in 
dynamically asymmetric mixtures show features reminis- 
cent of the cage effect and two-step relaxation (52[, one 
could imagine building two-step memory kernels with a 
plateau value around the above determined mn(q). 

How this should be done in practice will depend on 
the details of the dynamical theory at hand [56], [66| and 
is beyond the scope of the present paper. At this point, 
however, a few words of caution might be in order. First, 
the above argument specifically deals with one possible 
source of dynamical asymmetry, namely a wide separa- 
tion in the kinetic parameters (masses or free-diffusion 
coefficients) of the mixture components [47H49. 52]. It is 
not clear if and how it can be transferred to cases where 
the dynamical asymmetry has another origin, such as 
size disparity in dense systems (for recent studies, see 
Refs. [67|, |68[). Second, it is already known that not all 
theoretical frameworks are able to bridge between an- 
nealed mixtures and fluid-matrix systems. A prominent 
example is the mode-coupling theory [20H23| whose equa- 
tions for these two problems have been shown to be gener- 
ically incompatible [69]. In such a case, the connection 
described above does not bring any additional insight, 
neither into the dynamics of asymmetric bulk mixtures 
nor into that of RP and PP systems, and each problem 
has to be analyzed separately within the chosen theoret- 
ical framework. 

Obviously, and independently of the connection sug- 
gested above, the dynamics of the RP and PP systems 
has an interest of its own, as a simple model of a fluid 
in a disordered confining environment. In particular, be- 
cause the computation of static averages is made easy and 
well controlled in these systems thanks to the identities of 
Sec. IIIIl they are very appealing models to test dynamical 
theories which, quite gene rally and as illustrated by the 
mode-coupling theory [20H23I ] , require prior knowledge of 
a number of static properties to deliver their predictions. 
We stress here that the mode-coupling theory for fluids 
in random porous matrices [20- 23] does not make any as- 
sumption on the statistics of the disordered solid (beyond 
statistical homogeneity) and thus applies as it stands to 
RP and PP systems. 

These models can also be put to good use to investigate 
the configurational consequences of permanent trapping 
in systems with hard-core fluid-matrix interactions. In- 
deed, it is well known that, for such systems at any nonva- 
nishing matrix density, there is always a finite probability 
that particles are trapped in finite domains disconnected 
from the rest of the matrix void space [7(1 [Zl| • This trap- 
ping phenomenon culminates in a diffusion-localization 



transition which coincides with the percolation transition 
marking the point after which the matrix void space only 
consists of finite disconnected domains [71HZ3]. Thus, 
the exploration of the volume accessible to the fluid is 
never completely ergodic and one expects to see at least 
quantitative differences in configurational properties de- 
pending on whether the thermal averaging is performed 
as an ensemble or a trajectory average, denoted by (• • • ) 
and [•••], respectively. For a generic fluid-matrix system, 
such an investigation is hampered by a number of tech- 
nical diffi culties. For instance, if double averages of the 
form (A) and [A] are to be compared, the estimation of 
the former requires the use of smart simulation methods 
(typically, with nonlocal updates) in order to overcome 
the lack of er godicity of th e syste m. In the case of quan- 
tities such as (A)(B) and [A] [5], one has in addition to 
face the issues mentioned at the end of Sec. IIIIl for the 
computation of such correlation functions. Tackling the 
problem with a RP or a PP system does not eliminate 
all difficulties, but significant simplifications occur which 
should make the various computations much more acces- 
sible. Thus, thanks to Eq. (J7J), the double average (A) 
can be efficiently retrieved from studies of the bulk fluid 
system on which the fluid-matrix model is based, in which 
no trapping occurs, hence with no need for special sam- 
pling techniques. Also, rather than computing (A)(B) 
and [A][B], one can follow Eq. (f23|) and concentrate on 
correlation functions involving the positions of the fluid 
particles at the time of the preparation of the matrix, 
which are much easier to compute and have essentially 
the same physical content. 



VI. CONCLUSION 

In this paper, the statistical properties of the partly 
pinned fluid systems, a special class of models of fluids 
confined in disordered porous matrices, have been stud- 
ied. These systems, which in most respects are just typi- 
cal examples of fluid-matrix models, are singled out by a 
peculiar complementarity between the mobile and immo- 
bile components, which originates from the fact that the 
confining random solid is prepared in presence of and in 
equilibrium with the adsorbed fluid. A special symmetry 
results, with major repercussions on the configurational 
properties of the system. In particular, simple identities 
hold, which relate different types of configurational aver- 
ages, either relative to the fluid-matrix system or to the 
bulk fluid from which it is prepared. 

By taking advantage of these identities in computer 
simulation studies, it seems that, with these systems, in- 
teresting opportunities are opening up for efficient and 
accurate computations of quantities that are usually 
quite difficult to estimate in generic fluid-matrix mod- 
els, in particular the crucially important connected and 
disconnected correlation functions. In our opinion, this 
should make the PP fluid systems the models of choice for 
a number of theoretical developments in the field of fluids 
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adsorbed in disordered porous solids. Indeed, in many 
circumstances, either theories aim at predicting struc- 
tural properties, as do integral equation theories, and 
then it is clearly an advantage that their predictions can 
be readily tested against high-quality simulation data, 
or they take structural data as input in order to predict 
other physical properties, such as the dynamics, and then 
it is extremely convenient that this input can be easily 
and accurately determined. In all cases, one can reason- 
ably expect that the insight gained on the PP systems 
will be transferable to other types of fluid-matrix models, 
so that global progress can be made in our understanding 
of fluids in random environments. 

There is however a price to pay for these conveniences. 
Indeed, by construction of the model, one cannot change 
the parameters describing the adsorbed fluid without a 
similar and complementary modification of the matrix. 
This is clearly at odd with the typical experimental situ- 
ation, where, usually, the confining conditions are given 
and fixed and the fluid characteristics can be altered 
freely and independently, and might lead to difficulties 



of interpretation, when evolutions of the properties of 
the system across its parameter space are considered. 

Finally, some extensions of the present work could be of 
interest. One obvious direction is towards heterogeneous 
partly pinned fluid systems, in which all the particles of 
a bulk system located in a predefined spatial domain are 
stopped, leaving a fluid confined in a pore or in contact 
with a solid interface. This setup has been introduced 
in order to study various aspects of glass formation p78l — 
l8l| and an identity analogous to Eq. (TIT)) is known to 
hold (78|. Note however that, in order to prevent the in- 
vasion of the solid by the fluid, the authors of Ref. [78[ 
mention the need to introduce a hard separation at the 
fluid-matrix interface which might lead to complications 
and approximations absent in the homogeneous case. An- 
other possibility would be to translate the problem in the 
language of magnetism, where the pinning process would 
operate on spin variables. This would generate examples 
of spin systems with correlated random field and/or dilu- 
tion for which suitable adaptations of Eqs. (IT71) and ([25)1 
would apply. 



L. D. Gelb, K. E. Gubbins, R. Rad- 
hakrishnan, and M. Sliwinska-Bartkowiak, 

|Rep. Prog. Phys. 62, 1573 (1999) | 

M. Alcoutlabi and IT B. McKenna, 
|J Phys.: Condens. Matter 17, R461 (2005)| 
C. Alba-Simionesco, B. Coasne, G. Dosseh, G. Dudziak, 
K. E. Gubbins, R. Radhakrishnan, and M. Sliwinska- 
Bartkowiak, |J. Phys.: Condens. Matter 18, R15 (2006)| 
E. Kierl ik, M. L. Rosinberg, G . Tarjus, and P. Viot, 
|Phys. Chem. Chem. Phys. 3, 1201 (2001)| 
E. Kierlik, P. A. Monson, M. L. Rosinberg, L. Sarkisov, 
and G. Tarjus, |Phys. Rev. Lett. 87, 055701 (2001)| 
E. K ierlik, P. A. Monson, M. L. Rosinberg, and G. Tar- 
jus, [TTh^sT^ondensr^^teFli, 9295 (2002)] 



M. L. Rosinberg, E. Kierlik, 
|Europhys. Lett 62, 377 (2003) | 
R. L. C. Vink, K. Binder, 
|Phys. Rev. Lett. 97, 230603 (2006) | 



Vink, K. Binder, 



and G. Tarjus, 
and H. Lowen, 
and H. Lowen, 



R. L. 

|J. Phys.: Condens. Matter 20, 404222 (2008)] 

G. P ellicane, R. L. C. Vink, C. Caccamo, and H. Lowen, 

|J. Phys.: Condens. Matter 20, 115101 (2008)[ 

P. G. De Sanctis Lucentini and G. Pellicane, 



Phys. Rev. Lett. 101, 246101 (2008) 



P. Gallo, R. Pellarin, 
^Europhys. Lett. 57, 212 (2002)] 
P U^o, R Pellarin, 


and 
and 


M. 
M. 


Rovere, 
Rovere, 


[Phys. Rev. E 67, 041202 (2003) | 








P (JaUo, ~R. Pellarin, 


and 


M. 


Rovere, 


|Phys. Rev. E 68, 061209 (2003) 








A. Attili, P. Gallo, 


and 


M. 


Rovere, 


J. Chem. Phys. 123, 174510 (2005) 






P. Gallo, A. Attili, 


and 


M. 


Rovere, 



K. Kim, Europhys. Lett. 61, 790795 (2003) 
K. Kim, K. Miyazaki, and 



Saito, 



[20 
[21 

[22 
[23 
[24 

[25 

[26 

[27 

[28 

[29 
[30 

[31 

[32 

[33 
[34 
[35 
[36 

[37 

[38 
[39 

[40 

[41 

[42 

[43] 



Europhys. Lett. 88, 36002 (2009) 
V. Krakoviack, |Phys. Rev. Lett. 94, 065703 (2005) 
V. Krakoviack, |J. Phys.: Condens. Matter 17, S3565 (2005) 
V. Krakoviack, [ Phys. Rev. E~7 5, 031503 (2007) 
V. Krakoviack, Phys. Rev. E 79, 061501 (2009) 
J. Kurzidim, D. Coslovich, and G. Kahl, 
|Phys. Rev. Lett. 103, 138303 (2009)] 
W. G. Madden an d E. D. Glandt, 
|J. Stat. Phys. 51, 537 (1988)] 

L. A. Fanti, E. D. Glandt, and W. G. Madden, 
|J. Chem. Phys. 93, 5945 (1990)| 



W. G. Madden, J. Chem. Ph ys. 96, 5422 ( 1992) 
D. 



D. 



Glandt, 



M. Ford and E. 
|J. Chem. Phys. 100, 2391 (1994)| 
J. A. Given and 

|J. Chem. Phys. 97, 4573 (1992) 
E. Lomba, J. A. Given, G. Stell, J. J. Weis, and 



Stell, 



D. Levesque, [Phys. Rev. E 48, 233 (1993)] 

J. A. Given and G. Stell, |Physica A 209, 495 (1994)] 
M. L. Rosinberg, G. Tarjus, and G. Stell, 

| J. Chem. Phys. 100, 5172 (1994) | 

M. Schmidt, |Phys. Rev. E 66, 041108 (2002)| 

M. Schmidt, [j~Phys.: Condens. Matter 17, S3481 (2005) 
M. Schmidt, ]Pnys. Rev. E 79, 031405 (2009) | 
L. Lafuente and J. A~. Cuesta, 
|Phys. Rev. E 74, 041502 (2006) | 

P. R. Van Tassel, J. Talbot, P . Viot, and G. Tarjus, 

|Phys. Rev. E 56, R1299 (1997) | 

P. R. Van Tassel, |J. Chem. Phys. 107, 9530 (1997)1 
S. L. Zhao, W. Dong, and Q. H~~ 
| J. Chem. Phys. 127, 144701 (2007) | 
S. L. Zhao, W. Dong , and Q. H. 
|J. Mol. Liq. 136, 241 (2007)" 



Liu, 



Liu, 



P. R. Van Tassel, Phys. Rev. E 60, R25 (1999) 
L. Zhang and P. R. Van 
| J. Chem. Phys. 112, 3006 (2000) | 
L. Zhang and P. R. Van 



Tassel, 



Tassel, 



12 



son, |J. Chem. Phys. 106, 264 (1997) 



0>f(q,*))Pm(-q,0) 



(p f (q,0)>p m (-q,0) 



|Mol. Phy s. 98, 152 1 (2000)| 

[44] L. Zhang, S. Cheng, and P. R. Van Tassel, [64] This results from equalities such as (p f (q,i)p m (-q, 0)) 

|Phys. Rev. E 64, 042101 (2001) | 

W. G. Madden, |J. Chem. Phys. 102, 5572 (1995)1 
V. Krakoviack, E. Kierlik, M. L. Rosinberg, and G. Tar- 
jus, ^7^^7P^7ll5^TuW^^\ 
G. Viramontes-Gamboa, T. L~! Arauz-Lara, 



Arauz-Lara, 



and 

M. Medina-Noyola, [Phys. Rev. Lett. 75, 759 (1995) | 
G. Viramontes-Gamboa, M. Medina-Noyola, and J. L. 
Arauz-Lara, |Phys. Rev. E 52, 4035 (1995) | 
M. A. Chavez-Rojo, R. Juarez-Maldonado, and 
M. Medina-Noyola, Phys. Rev. E 77, 040401 (2008) | 
R. Chang, K. Jagannathan, and A. Yethiraj, 
|Phys. Rev. E 69, 051101 (2004)| 

J. Mittal, J. R. Errington, and T. M. Truskett, 
|Phys. Rev. E 74, 040102 (2006)| 

W. Fenz, I. M. Mryglod, O. P rytula, and R. Folk, 
|Phys. Rev. E 80, 021202 (2009)1 

G. Cruz de Leon, J. M. Saucedo-Solorio, and J. L. Arauz- 
Lara, |Phys. Rev. Lett. 81, 1122 (1998)| 
G. Cruz de Leon and T. L. 
|Phys. Rev. E 59, 4203 (1999) | 

Multicomponent versions are possible for which every 
species has its own pinning fraction. The generalization 
of the results of the present work to this case is straight- 
forward. 

J. -P. Hansen and I. R. McDonald, Theory of simple liq- 
uids, Third edition (Academic Press, London, 2006). 
I. M. Lifshits, S. A. Gredeskul, and L. A. Pastur, Intro- 
duction to the theory of disordered systems (Wiley, New 
York, 1988). 

A similar result for functions of the fluid variables only 
has been reported for heterogeneous partly pinned sys- 
tems in Ref. [7^. 

A. Meroni, D. Levesque, and J. J. Weis, 
|J. Chem. Phys. 105, 1101 (1996)| 

D. F. Schwanzer, D. Coslovich, J. Kurzidim, and 
G. Kahl, |Mol. Phys. 107, 433 (2009)| 

C. 



G. I. Menon and 
|Phys. Rev. Lett. 73, 1023~( 1994) 



Dasgupta, 



C. Caccamo, |Phys. R ep. 274, 1 (1996) 
E. Kierlik, M. L. Rosinberg, G. Tarjus, and P. A. Mon- 



[65 
[66 
[67 
[68 
[69 

[70 
[71 
[72 
[73 
[74 

[75; 

[76 
[77 
[78 
[79 
[80 
[81 



{p { (q, 0)p m ( 
G. N agele, 

j 



0)). 

Bergenholtz, 



and J. K. G. Dhont, 



Chem. Phys. 110, 7037 (1999) 
P. Boon and S 



New York, 
A. J. 



Yip, Molecular hydrodynamics (Dover, 
1991). 

Moreno 



and 



Colmenero, 



|Phys. Rev. E 74, 21409 (2006) 

T. Voigtmann and J. Horbach, 

|Phys. Rev. Lett.103, 205901 (2009)] 
Following Ref. [H, the likely reason for this is that the 
derivation of the mode-coupling theory for fluid-matrix 
systems critically relies on the distinction between con- 
nected and blocking correlation functions that does not 
make sense for bulk systems. 

J. Kertesz and J. Metzger, J. Phys. A: Math. Gen. 16, 
L735 (1983). 

A. Kammerer, F. Honing, and T. Franosch, Europhys. 
Lett. 84, 66002 (2008). 

F. Honing, T. Franosch, and E. Frey, 

|Phys. Rev. Lett. 96, 165901 (2006) | 
F. Honing, T. Munk, E. Frey, and T. Franosch, 
|J. Chem. Phys. 128, 164517 (2008) | 

B. J. Sung and 
|Phys Rev. Lett. 96, 228T03~( 2006) 
B. J. Sung and A. Yethiraj, J. Phys. Chem. B 112, 143 
(2008). 

B. J. Sung and 

|J. Chem. Phys .128, 054702~( 2008) 
S. Babu, J. C. Gimel, and T. Nicolai, J. Phys. Chem. B 
112, 743 (2008). 

P. Scheidler, W. Kob, and 
|J. Phys. Chem. B 108, 6673 (2004)"] 
J.-P. 

|J. Chem. Phyl 

A. Cavagna, T. ST Grigera, and P. 



A. 



A. 



Bouchaud 



and 



121, 7347 (2004) 



Yethiraj , 



Yethiraj , 



K. Binder, 



Biroli, 



Verrocchio, 



|Phys. Rev. Lett. 98, 18780T7 2007) 

G. Biroli, J.-P. Bouchaud, A. Cavagna, T. S. Grigera, 
and P. Verrocchio, |Nature Phys. 4, 771 (20"08)| 



